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■ Abstract 

, We have calculated the form-factors F and G in K — > tyty£u decays (Ku) to two-loop 

order in Chiral Perturbation Theory (ChPT). Combining this together with earlier two- loop 

ON ', calculations an updated set of values for the L\, the ChPT constants at 0(p ), is obtained. 

CN ■ We discuss the uncertainties in the determination and the changes compared to previous 

estimates. 

(N ' 
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OO ! 

, 1. The theoretical study of A — > -khIv decays, Ku decays, provides an interesting possibility to 

C*~) ' test our understanding of the long-distance dynamics of the strong sector in the Standard Model. 

CN ! Perturbative QCD can not be applied at energies lower than the spontaneous symmetry scale, 

A x ~ m p . One has to resort to an effective approach to the full theory, Chiral Perturbation 
Theory (ChPT) 0, to obtain a reliable model independent description of physical processes. The 
main purpose in the study of Ku decays is twofold: i) This decay is one of the cleanest sources 
of 7T7T pairs at low-energies and thus provides us with the possibility to check 7T7r-scattering near 
threshold, ii) The form-factors themselves are also directly of interest. They provide a direct test 
Cln. of our understanding of the three-flavour sector and are in addition one of the major inputs to 

determine the needed constants to predict other quantities. The early history can be found in the 
review M. Ku have been treated in the context of ChPT at one-loop || and in a dispersively 
improved one- loop representation [Q . These studies delivered the standard values of the low-energy 
\ constants L\,I7^ and L^. 

5_i ■ On the other hand, 7T7r-scattering has been studied at two-loop order in two-flavour ChPT, 

direct Q and via dispersive methods The main remaining uncertainty on the prediction of 
7T7r-scattcring is the value of the low-energy constants. In Q the standard values from [Q were 
used. A Roy equation analysis using high energy 7T7r-data, gave significantly different values for 
these constants [frl. There are several possible sources for this discrepancy: the Omnes representa- 
tion of Q is not sufficient for Ki±; the presence of large systematic effects in the 7T7r-scattering data 
base and Roy analysis; or large corrections to the one-loop relation between two- and three-flavour 
ChPT constants. Therefore, and since new experiments are underway, a full two-loop calculation 
of Ku is necessary. A first step was the estimate of the double logarithm corrections of 0(p 6 ) to 
and Ku ||. It was shown there they could be large. In addition pions in the 1 = 0, S-wave 
have strong final state interactions. We describe here first results of the full Ku C(p 6 ) calculation. 

2. Ku decays are described by the matrix element 

M ij = ^V: s tt(p„) 7/ ,(l - l 5 )v(Pi)(7r' l ( Pl )7rHp 3 )\V^ - A»\K(k)) . (1) 

Lorentz invariancc allows to parametrize the hadronic part by four form-factors 

(ttVI^IA) = — \P^F^ + Q^G lJ + L^Rv] , 

m K 

(7rV|V|tf) = -^e^L u P a Q^ , (2) 



K 



1 




Figure 1: Vertex topology. Dots refer to strong vertices or current insertions. 



with P = pi + pj , Q = pi — pj and L = pi + p v . Pi is the momentum of 7r\ All four form- factors 
are dimensionless functions of three variables, s n = P ■ P,u n = (k — pi) 2 and t n = (k — Pj) 2 - 
They can be decomposed, into a part symmetric or antisymmetric under t n <-> interchange, 
corresponding respectively to the isospin or 1 part for F, R and isospin 1 or for G, H . These 
processes are dominated by F and G. H has a fairly small influence and a one-loop calculation, as 
performed in 0, is sufficient to reach 0(p e ) since it comes through the odd-intrinsic-parity sector. 
R only contributes proportional to m 2 and can thus be neglected in the decays with an electron. 

In the remainder we concentrate on the ij = H — channel. The others can be derived from it 
using isospin relations. The older available experiments are compatible with the K + — > n + Tr~ e + v e 



experiment |10|. The most recent other measurement is Kj, — > Tv^vPe^v 

The experiment relied on a partial wave analysis of the form-factors. The analysis only 
kept s and p waves, the effect of (i-waves as well as of si was within the measurement errors and 
neglected in the further analysis. They used the parametrization 

F = f s e iS ° + f p e tS " cos W + D-wavc + . . . , G = ge lS " + D-wave + . . . , (3) 

where Si are the i-wave strong two-pion final state phase shifts, and / s , f p and g are defined to 
be real functions with argument s n . They then divided the data in 5 energy bins, observed an 
dependence in f s and one compatible with it for g. The final fit was performed using 

f s (q 2 ) = f s (0) (1 + \ f q 2 ) , g(q 2 ) = g(0) (l + \ g q 2 ) , q 2 = (s„ - 4m*)/(4m») . (4) 

with Xj = X g . Furthermore, the f p value was compatible with zero. The results obtained, with 
sin 6c = 0.22, are 

/ s (0) = 5.59 ±0.14, 9 (0) = 4.77 ±0.27, X f = 0.08 ± 0.02 . (5) 

Using isospin the result of jllj is g(0) = 5.50 ± 0.50 which using PDG procedures leads to a 
combined 

g(0) = 4.93 ± 0.31 . (6) 



Note that |11 neglected f p which can account for up to 5% of the value of g(0). 



3. As mentioned above, earlier calculations gave an indication of potentially sizeable corrections 
at 0(p 6 ). In order to obtain the full correction we have evaluated both form factors, F and G, 
at next-to-next-to-leading order (NNLO) in ChPT following a diagrammatic approach. One faces 
the evaluation of a large number of diagrams but only one new topology, shown in Fig. |l], appears 
w.r.t. the vector- vector and axial- vector-axial- vector two point-functions fl3] ]. It involves a new 
set of integrals, the vertex-integrals, that can be obtained in terms of a two-parameter integral 
representation |Q. 

The full evaluation leads to rather long expressions. More than half of the complexity stems 
from the topology of Fig. [ljdue to the large number of possible mass and momentum combinations. 
We have performed several checks on our full expressions: i) As a basic test of the algebraic 
programs we recover the one- loop expressions |3|, ^j. ii) All non-local divergences cancel when 
adding the full set of diagrams together with wave function renormalization. Hi) The polynomial 
divergences also cancel against the counter-terms determined in general for the even-intrinsic parity 
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C(p 6 ) sector iv) The double logarithms terms are in agreement with those of M. v) Some 
diagrams involving three-point one-loop integrals can also be obtained via the renormalization 
of the one-loop graphs, and hence the three-point integrals should cancel in the final result, vi) 
The sunset-type integrals encountered in Jl4[ are in agreement with those already calculated in 
p3[ using different methods, vii) The two-particle discontinuities of the vertex-type integrals have 
been checked using the Cutkosky rules below the three-particle thresholds, viii) Our results satisfy 
the isospin relation M+~ = M~%/2 + M oa . In view of these checks, we trust our calculations of 
the matrix-elements. 

To obtain the final answer, we shift the bare quantities -masses and decay constants- to the 
renormalized physical ones and add the part coming from the 0(p e ) Lagrangian as determined in 
& 

4. One of the main problems is how to deal with the coefficients of the 0(p 6 ) Lagrangian. 
Here we follow the philosophy of most of the other 0(p 6 ) calculations and estimate them using 
resonance exchange. This worked well at 0(p 4 ) fl7fl . In the second paper of jL8| this was found 
to agree reasonably well at 0(p 6 ) for the form- factors considered there. We use the notation as 
defined in fl3|| but include also terms which did not contribute there. Specifically we use for the 
vector nonet matrix 

+ia v (V li [u v ,f fW ]) (7) 
and for the axial-vector nonet 

£a = -\{A, V A^) + \m 2 A {A^) - -^{A^D+l^ (A^u")Wa ) (W> ■ (») 

We have chosen a representation where the vectors and the axial-vectors only start contributing 
to the mesonic Lagrangian at 0(p e ). For the scalar nonet we take 

£s = ~<V M SV"S) - \m 2 s {S 2 ) + c d (Su^) + c m (S X +) ■ 

We fixed the parameters as much as possible from experiment or the comparison with 0(p i ) j^, [l7| . 
For the remainder we use the values predicted by the Nambu-Jona-Lasinio model ]isj . The spe- 
cific inputs used are my = 0.77 GeV, m A =1.23 GeV, m s =0.98 GeV, fy = 0.20, f x = -0.025, 
g v = 0.09, ay = -0.014, f A = 0.1, 7^ = 0.006, 1a 2) = -0.01, c m = 42 MeV and c d = 32 MeV. 



These are used in all the fits of Table |1], except fit 6 where we kept only fy , f x and gy as nonzero 
couplings. 

5. ChPT in the meson sector has as parameters at 0(p 2 ) F , B and the quark masses rh = 
m u — rrid and m s . The latter only appear multiplied with Bq for a total of 3 parameters. At 0(p 4 ), 
we have L\, . . . , L\ Q Q and at 0(p 6 ) there are 90 additional parameters [[l6]]. In the two-flavour 
sector the equivalent of Lg and L\ Q has been determined to NNLO [|l8). L\ has no influence on 
the other quantities considered here and the influence of Lg is small. We will therefore use the 
standard values of [@, [§. For the 0(p 6 ) we use the resonance saturation approximation described 
above. 

As input we use = 92 A MeV, Fk/F^ — 1.22 ± 0.01 and the physical masses of pion, 
kaon and eta |Q in addition to / s (0), 5(0) and A/. The final input is the ratio m s /rh. We 
further assume that the parameters L\ and Lg are small because they vanish in the large 7V c -limit. 
When calculating the form-factors for K + — > 7r + 7r~e + t/, F^ and Fk we use the physical charged 
masses m T + = 139.56995 MeV; rriK+ = 493.677 MeV and m v = 547.30 MeV. When comparing 
with the quark-mass ratios we use the 7r°-mass, m v a — 134.9764 MeV, the physical 77-mass and 
m 2 K = l/2(m 2 <+ + m 2 K0 — 1.8(m 2 + — m 2 )) = (494.53 MeV) 2 . This removes the electromagnetic 
mass corrections including estimated corrections to Dashen's theorem [BM. 
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Table 1: Results for L^(^l) for the various fits described in the main text. Notice that LJ, Lg and 
Lg are input, L\ = 0, Lg — and Lg = 6.90 10~ 3 . Errors are fitting errors described in the text. 
All L£(/j) values quoted have been brought to the scale /i = 0.77 GeV. The standard values are 
m s /m = 24, ^ = 0.336 GeV and s e = 0. 





Main Fit 


[1, 4| 


0(p 4 ) 


fit 2 


fit 3 


fit 4 


fit 5 


fit 6 


fit 7 


fit 8 


fit 9 


10 3 L{ 


0.52±0.23 


0.37±0.23 


0.46 


0.52 


0.50 


0.49 


0.52 


0.45 


0.42 


0.63 


0.65 


10 3 LI 


0.72±0.24 


1.35±0.23 


1.49 


0.73 


0.67 


0.74 


0.80 


0.51 


1.05 


0.73 


0.85 


10 3 LI 


-2.69±0.99 


-3.5 ±0.85 


-3.18 


-2.70 


-2.57 


-2.73 


-2.74 


-2.15 


-2.91 


-2.68 


-3.27 


10 3 U h 


0.65±0.12 


1.4 ±0.5 


1.46 


0.63 


0.65 


0.65 


0.82 


0.67 


0.99 


0.51 


0.60 


10 3 U, 


-0.31±0.15 


-0.4 ±0.2 


-0.49 


-0.24 


-0.31 


-0.31 


-0.35 


-0.30 


-0.28 


-0.29 


-0.31 


10 3 LI 


0.48±0.18 


0.9 ±0.3 


1.08 


0.36 


0.48 


0.49 


0.59 


0.48 


0.52 


0.45 


0.49 


changed 






0(p 4 ) 


m a /rh 


s/si 




rr.rr 

L14, -kg 


V M 






5(0) 


quantity 








26 


0.1 


0.293 


-0.3; -0.2 


only 


0.5 


1.0 


4.93 


Unit 










GeV 


GeV 


10" 3 




GeV 


GeV 





In order to compare with expression (|4j) we have to pull out some of the kinematic behaviour 
encoded in the matrix elements. We set 

F(s w ,s e , cos On = 0) , 3(0) = G(s n ,Si, cos W = 0) , 
( \F( 3 ' v , Se ,caa0„=O)\ \ 4m 2 (Q) 
\\F(s„,st,co80„=0)\ J s'„-s w ' [ ' 

We use s-n = (2771^ ± 1 MeV) 2 to avoid numerical integration problems very near threshold, 
s'^ = (336 MeV) 2 and sg = 0. / s (0), g(0) and A/ are fitted with the experimental errors quoted in 

®- ™ 

Fk/F^ is calculated with the formulas of |13J and required to be 1.22 ± 0.01. The lowest order 

meson masses, m^, are calculated from the physical ones using the formulas of |l3[ ] and we require 
m s /m = {2ml K - m^/m 2 ^ = (3m^ - m^/Qm^) = 24 with an error of 10%. 

The errors quoted above are the source of the errors quoted in our main fit in Table @. They do 
not include an estimate of the theoretical uncertainties due to higher orders. The error procedure 
is the same as in j|, since we use the same data the errors have thus not changed. 

We first display our main fit and the canonical values. The 0(p 4 ) fit illustrates the change 
due to the 0(p 6 ) corrections. It is similar to the one-loop fits performed in the difference due 
mainly to using L\ — and Lg = 0. The remainder of the table shows the variation with the 
choices of input. Fit 2 is with a different choice of input quark-mass ratio, m s /m — 26, fit 3 with 
a different choice of y/sj = 0.1 GeV, fit 4 of 4, fit 5 with L\ = -0.3 lO" 3 and L r 6 = -0.2 10~ 3 . 
Setting LI — does not change results within the accuracy displayed. Fit 6 is with only the 
vector meson part in the resonance estimate. As can be seen the main contribution comes from 
the vector-exchange with the experimentally determined quantities. The rest is from the scalar 
sector. The axial-vectors contribute only for sg 7^ 0. Fit 7 and 8 are with two different values of 
the renormalization scale, /i = 0.5 GeV and /j, = 1.0 GeV. Our full result is /i- independent as it 
should be, the difference is due to the fact that we estimated the value of the 0(p 6 ) coefficients 
at the scale \x. In fit 9 we present the result with the g(0) value of |ll| included. 

6. In conclusion, we have performed a NNLO calculation of the F and G form-factors in 
Kg4 decays. We then used this calculation together with earlier calculations of masses and decay 
constants to update the ChPT parameters at 0(p 4 ) with a certain range of choices for the fitting of 
the available K n data and the estimates of the 0(p 6 ) parameters. None of the variations of input 
change the parameters outside the experimentally determined error except the rather extreme case 
of p — 0.5 GeV. 

When looking at the calculation it can be seen that for the form-factor F the estimate 
of higher orders of 0] was in the right direction. The changes are somewhat larger than naively 
expected. 



fs(0) = 
A/ = 
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L\ versus L 3 

1.1 I 1 1 1 1 1 r 
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L' 2 versus L 3 

1.4 | , , 1 , , r 




1 1 1 1 1 1 1 1 1 1 1 

-5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 
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Figure 2: Correlation plot of L\ versus L 3 and L r 2 versus L 3 . Shown are the points from a 
distribution of L\ with the \ 2 calculated from the observables described in the text as for our 
main fit. Only points within a 68% confidence level are included. 



The errors should be taken with caution, Uj and Lg are strongly anti-correlated and LI is 
strongly anti-correlated with L\ and I7 2 . As an example of this we have shown in Fig. |2j a 
distribution of sets of L\ that fall within the 68% confidence level limit. If we choose as fit 
variables X\ = L\- 2L\ - L r 3 , X 2 = L\ and X 3 = (L 2 - 2L[)/L r 3 (§, X 3 is little correlated with 
the others. We obtain 

X 3 = 0.121°;°?. (10) 

The large iV c -prediction |_X"3 ] <C 1 is obviously well satisfied. These questions will be discussed in 
more detail in a future publication. 

With these values for the low-energy constants we can see how the various quantities behave. 
The numbers correspond to 0(p 2 ), 0(p A ) and 0(p 6 ) contributions. 

m^o/m^o|exp = 0.740 + 0.007 + 0.253, 

m 2 K /m 2 K \exp = 0.689 + 0.024 + 0.287, 

m^/m2| ex p = 0.735- 0.066 + 0.331, 

F K /F„ = 1.000 + 0.135 + 0.085, 

Fn/F = 1.000 + 0.136 - 0.076, 

/ s (0) = 3.78+ 1.15 + 0.66, 

g(Q) = 3.78 + 0.82 + 0.17, 

A/ = 0.000 + 0.127- 0.047. (11) 

These numbers were calculated before rounding the L\ to the number of significant digits given 
in Table |l|. Contrary to what was observed in ^ ^l) the slope is now mainly from 0(p 4 )-effects. 

Making use of the central values of the main fit in Table |l| together with the 0(p 4 ) relations 
between the two- and three-flavour low-energy constants we estimate l\ and l 2 

h = 0.4(-0.2) , l 2 = 4.9(5.2) . (12) 

The values in brackets are those using the values of fit 9. We remind the reader that the relation 
used has possibly large 0(p 6 ) corrections. These should be compared respectively with —1.7 and 
6.1 using the L\ from {§ and with l 2 = 4.2 0. 

In conclusion, we have performed a full 0(p 6 ) calculation of Ku decays. We then performed a 
full refit of the 0(p 4 ) parameters of ChPT and discussed the changes and the validity of a large- N c 
prediction. 
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